Quadrupole anisotropy in dihadron azimuthal correlations in central d+Au collisions 

at ys^=200 GeV 
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The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) reports measurements 
of azimuthal dihadron correlations near midrapidity in d+Au collisions at ys„„=200 GeV. These 
measurements complement recent analyses by experiments at the Large Hadron Collider (LHC) 
involving central p+Pb collisions at ^s„j^=5.02 TeV, which have indicated strong anisotropic long- 
range correlations in angular distributions of hadron pairs. The origin of these anisotropies is 
currently unknown. Various competing explanations include parton saturation and hydrodynamic 
flow. We observe qualitatively similar anisotropies at RHIC to those seen at the LHC, and when 
both are divided by an estimate of the initial-state eccentricity, the anisotropies follow a common 
multiplicity scaling. This scaling is also found to extend to heavy ion data at RHIC and the LHC, 
where the anisotropies are widely thought to be due to hydrodynamic flow. The results presented 
here, at much lower collision energy and with a deuteron projectile (instead of a proton), provide 
important new information for understanding the origin of these new long-range correlations. 

PACS numbers: 25.75.Dw 



Proton- and deuteron-nucleus collisions at relativistic 
energies have been studied in order to provide baseline 
measurements for heavy ion collisions. In p((i)-f A col- 
lisions, initial-state nuclear effects are present; however, 
the formation of hot quark-gluon matter as created in 
heavy ion collisions is not commonly expected. Recently 
there has been significant interest in the physics of very 
high multiplicity events in small collision systems, mo- 
tivated by the observation of a small azimuthal angle 
(A0) large pseudorapidity (Ary) correlation of primarily 
low pt particles in very high multiplicity p-\-p collisions at 
7 TeV [l| . The correlation strikingly resembles the "near- 
side ridge" observed in Au-^Au at 200 GeV 0,13]. The 
initial_p-|-p result sparked considerable theoretical inter- 
est \^]^ . Recently, a similar effect was observed in p+Fh 
collisions at y/s^ = 5.02 TeV 0. Subsequent work re- 
moved centrality independent correlations (largely from 
jet fragmentation) by looking at the difference in correla- 
tions between central and peripheral events and has addi- 
tionally uncovered similar long-range A77 correlations on 
the opposite side (A^ w tt) beyond those expected from 
fragmentation of recoiling jets measured by ALICE Q 
and ATLAS [9|]. The effect appears as a longitudinally- 
extended azimuthal modulation with a predominantly 
quadrupole component (i.e. cos 2A(/)) and bears a quali- 
tative resemblance in both magnitude and pt dependence 
to elliptic flow measurements in heavy ion collisions, 
where the large quadrupole modulation is understood to 
be caused by the initial-state spatial anisotropy followed 
by a nearly inviscid hydrodynamic expansion [10| . A vari- 
ety of physical mechanisms have been invoked to explain 
the observed anisotropies in p-\-Ph including gluon satu- 



ration [gI- HiUi^ . hydrodynamics [1, [l3,[l3], multiparton 
interactions [16j, and final-state expansion effects 17 1. 

Previous analyses involving two-particle correlations 
from d-\-Au collisions at RHIC have not indicated any 
long-range features at small Acj) 0, [l8l - Eo| . However, 
these measurements involved pt selections which em- 
phasize jet-like correlations, rather than the underlying 
event. Also, Refs. 1^, 2^ were based on d-|-Au collisions 
recorded in 2003 with a small data sample which limited 
the statistical significance of the results. 

We present here the first analysis of very central d-l-Au 
events to measure hadron correlations between midrapid- 
ity particles at ^/s^ = 200 GeV. The center of mass 
energy per nucleon is a factor of 25 lower than at the 
LHC, and another potentially key difference is the use of 
a deuteron as the projectile nucleus rather than a pro- 
ton. In Ref. within the context of a Monte Carlo 
Glauber model, the calculated initial spatial eccentric- 
ity of the participating nucleons, £2, for central (large 
number of participants) d-\-Ph is more than a factor of 
two larger than in central p-|-Pb collisions at LHC en- 
ergies. We find the initial spatial eccentricity £2 from 
Monte Carlo Glauber 2l| for d-|-Au at RHIC energies to 
be similar to the d-l-Pb calculations at LHC energies. 

The results presented here are based on 1.56 billion 
minimum-bias d-|-Au collisions at ^s„„ = 200 GeV 
recorded with the PHENIX detector in 2008. The event 
centrality in d-l-Au is determined by categorizing the in- 
tegrated charge by upper percentile as seen by a beam- 
beam counter (BBC) facing the incoming Au nucleus [H]. 
Here we isolate a more central sample than previously 
analyzed, in order to compare more closely to the LHC 
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results. We use central and peripheral event samples 
comprising the top 5% and 50-88% of the total charge 
distributions, respectively. 

This analysis considers charged hadrons measured 
within the two PHENIX central arm spectrometers. 
Each arm covers nominally 7r/2 in azimuth and has a 
pseudorapidity acceptance of \ri\ < 0.35. Charged tracks 
are reconstructed using the drift chambers with a hit 
association requirement in two layers of multiwire pro- 
portional chambers with pad readout, achieving a mo- 
mentum resolution of 0.7% © l.l%p (GeV/c). Only 
tracks with full and unambiguous drift-chamber and pad- 
chamber- 1 hit information are used. Electrons are re- 
jected with a veto in the ring-imaging Cerenkov (RICH) 
counters. 

All pairs satisfying the tracking cuts within an event 
are measured. The measured pairs are then corrected 
for the PHENIX azimuthal acceptance through use of 
mixed event distributions. The conditional yield of pairs 
is determined by: 
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where iV* is the number of trigger hadrons (trigger 
hadrons are those which have the momenta required to 
begin the search for a pair of hadrons) and -/V^p^'J (-^mix ^) 
is the number of pairs from the same (mixed) events. 
Mixed pairs are constructed with particles from differ- 
ent events within the same 5% centrality class and with 
event vertices within 5 cm of each other. Since the fo- 
cus of this analysis is on the shape of the distributions, 
no correction is applied for the track reconstruction ef- 
ficiency, which has a negligible dependence on centrality 
for c?-|-Au track multiplicities. 

In order to make direct comparisons between our mea- 
surements and recent ATLAS p-t-Pb results 0, we fol- 
low a similar analysis procedure. Charged hadron selec- 
tions are made at different momenta from 0.5 through 
3.5 GeV/c. For this analysis, each pair includes at least 
one particle at low pT (0.5 < pr < 0.75 GeV/c) in or- 
der to enhance the sensitivity to the nonjet phenomena. 
The pairs are restricted to pseudorapidity separations 
of 0.48 < IA77I < 0.7, in order to minimize the con- 
tribution from small-angle correlations arising from res- 
onances, Bose-Einstein correlations, and jet fragmenta- 
tion. This pseudorapidity gap is chosen to be as large as 
possible within the PHENIX tracking acceptance, while 
still preserving an adequate statistical sample size. 

The associated yield due to azimuthally uncorrelated 
background is estimated by means of the zero-yield-at- 
minimum (ZYAM) procedure 23[. This background con- 
tribution is obtained for both the central and peripheral 
samples by performing fits to the conditional yields using 
a functional form composed of a constant pedestal and 
two Gaussian peaks, centered at A<j> = and tt. The 
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FIG. 1: (color online) Azimuthal conditional yields, Y ( A(j)), 
for (open [black] squares) 0-5% most central and (open 
[black] circles) peripheral (50-88% least central) collisions 
with a minimum Arj separation of 0.48 units, (filled [blue] 
circles) Difference AY{A(j}), which is ([blue] curve) fit to 
ao + 2a2 cos(2A(^), where ao and a2 are computed directly 
from the data, (shaded [blue] band) Statistical uncertainty 
on 02. No correction for the A(/!> independent reconstruction 
efficiency has been applied. 



minimum of this function, 6zyaMj is subtracted from the 
Acf) distributions, and the result is y(A(/)): 
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The conditional yields Yc{A(j)) and Yp{A(j)) (central and 
peripheral events, respectively) are shown in Fig.[Tl along 
with their difference AF(A0) = Yc{A(j)) - Yp{A(l}). As 
in Ref. [91], this subtraction removes any centrality inde- 
pendent correlations, such as effects from unmodified jet 
fragmentation, resonances and HBT. In the absence of 
any centrality dependence, Yc{A(j)) and Yp{A(j)) should 
be identical. Due to the limitations of our method, any 
signal in the peripheral events is subtracted from the cen- 
tral events. We see that Yc{A(j>) is significantly larger 
than Yp{A(t)) for Acj) near and tt. 

In a manner similar to Ref. 9] , we find that the differ- 
ence with centrality is well described by the symmetric 
form: AY{A<j)) w oq + 2a2 cos (2A0) as demonstrated 
in Fig. [TJ The coefficients a„ and their statistical uncer- 
tainties are computed from the AY{A(f>) distributions as: 
a„ = (Ay(A0)cos(nA0)). The cos(2A</)) modulation 
appears as the dominant component of the anisotropy 
for all trigger/partner combinations as will be quantified 
below. 

The PHENIX central arm spectrometers lack sufficient 
\Ari\ acceptance to completely exclude the near-side jet 
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FIG. 2: (color online) Sample comparison of Y (Acp) and 
Ay {A<j>) for same and oppositely charged pairs for 1.25 < 
Pt < 1.5 GeV/c and 0.48 < |A?7| < 0.7. The symbols, curve, 
and shaded band are as described in the Fig. [1] caption. 
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peak. To assess the systematic influence of any residual 
unmodified jet correlations, we analyzed charge-selected 
correlations. Charge-ordering is a known feature of jet 
fragmentation which leads to enhancement of the jet cor- 
relation in opposite-sign pairs, and suppression in like- 
sign pairs, in the near side peak (e.g. Ref. [2^). A rep- 
resentative Pt selection of ^(A^) and AY{A(j)) distribu- 
tions are shown in Fig. [21 where all charge combinations 
exhibit a significant cos 2 modulation. The magnitude 
of the modulation is larger in the opposite-sign case, indi- 
cating some residual unmodified jet correlation contribu- 
tion. We also varied the {Arjl window which changes the 
residual jet contribution. Both of these cross-checks are 
used to estimate the systematic uncertainty, as discussed 
later. 

In order to quantify the relative amplitude of the az- 
imuthal modulation we define c„ = a„/ (^Iyam + 0^0) 
where h 



ZYAM 



is 6zYAM in central events. This quantity is 
shown as a function of associated pr in Fig. |3] for central 
(0-5%) collisions. 

The centrality dependence will be analyzed in further 
detail in a forthcoming publication, though we note that 
we have observed a signal of similar magnitude for the 
0-20% most central colhsions. The ATLAS C2 results [9] 
have a qualitatively similar dependence, but with a 
significantly smaller magnitude. However, it must be 
noted that the C2 values from PHENIX and ATLAS are 
not directly comparable since C2 is a function of the pT 
of both particles and the trigger particle pT range is not 
identical in the two analyses. ATLAS has also used a 
much larger Arj separation between the particles. 

The C3 values, shown in Fig.|3l are small relative to C2. 
Fitting the C3 data to a constant yields (6±4)x 10^"* with 
a x^/dof of 8.4/7 (statistical uncertainties only). The 
current precision is inadequate to reveal the existence of 
a significant C3 signal. 

In p+Ph collisions the signal is seen in long range A77 



FIG. 3: (color online) The nth-order pair anisotropy, c„, of 
the central collision excess as a function of associated par- 
ticle p^. PHENIX (fiUed [red] circles) C2 and (open [black] 
circles) C3 are for 0.5 < Pt < 0.75 GeV/c, 0.48 < |A»7| < 0.7 
and ATLAS (filled [green] squares) C2 [9[ are for 0.5 < Px < 
4.0 GeV/c, 2 < \Ar]\ < 5. 
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FIG. 4: (color online) Charged hadron second-order 
anisotropy, «2, as a function transverse momentum for (filled 
[blue] chcles) PHENIX and (open [black] circles) ATLAS ^. 
Also shown are a hydrodynamic calculation [13, [1^ for (up- 
per [blue] curve) d-l-Au collisions at y^s]^ = 200 GeV 
and (lower [black] curve) 0-4% central p-|-Pb collisions at 



= 4.4 TeV. 



correlations. Here, signal is measured at midrapidity, 
but it is natural to ask if previous PHENIX rapidity sep- 
arated correlation measurements [l8| would have been 
sensitive to a signal of this magnitude, if it is present. 
The maximum C2 observed here is approximately a 1% 
modulation about the background level. Overlaying a 
modulation of this size on the conditional yields shown in 
Fig. 1 of Ref. 18| shows that the modulation on the near 
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side is small compared with the statistical uncertainties 
on the points. In the current analysis, both particles are 
near midrapidity, while the analysis in Ref [18| includes 
one of the particles very forward (3.0 < rj < 3.8) in the 
d-going direction. Thus, with the current results we can- 
not determine whether the signal observed here persists 
for r/ >3. 

A measure of the single-particle anisotropy, V2 , can be 
obtained under the assumption of factorization (i^-Sl, 
which gives the relation C2 (pt^Pt) — ^2 {Pt) ^ "2 {Pt)- 
We have varied ptp and recomputed V2 (pt) and find no 
significant deviation from this factorization hypothesis. 
The calculated single particle V2 is shown in Fig. |4l and 
also compared with the ATLAS [§] results, again reveal- 
ing qualitatively similar pT dependence with a signifi- 
cantly larger magnitude. We also compare the V2 re- 



Q g ^ cm PHENIX d+Au ((s=200 GeV) 



suits to a hydrodynamic calculation [ij, |25| and find 
good agreement between the data and the calculation, 
which predicts larger anisotropy in d-|-Au than p-\-Ph col- 
lisions (the calculation for p-|-Pb is for 0-4% centrality at 
4.4 TeV, not 0-2% central at 5.02 TeV as in the data). 

The systematic uncertainties as shown in Figs. [3] and 
m are estimated as the root-mean-squared variation of 
the same-sign and opposite-sign c„ measurements about 
the combined value to reflect the influence of possible re- 
maining jet correlations. This systematic uncertainty is 
applied symmetrically, since the influence of the jet con- 
tribution is not known. As a test, the A77 interval was 
varied from the nominal value of 0.48 to 0.36 and 0.60. 
The Cn values remained unchanged within statistical un- 
certainties, with the qualiflcation that the jAyyl > 0.6 
sample lacks sufficient statistics for a precise comparison 
at higher px- We also produced V2 values with different 
trigger particle momentum selections and found no sig- 
nificant change in the extracted values. Other sources 
of uncertainty, such as occupancy and acceptance correc- 
tions, were also found to have negligible effect on these 
results. 

In order to further investigate the origin of this effect 
in Fig. [5] we plot the RHIC and LHC results scaled by 
£2 as calculated in a Glauber Monte Carlo as a func- 
tion of the charged particle multiplicity at midrapidity. 
The 0-5% d-l-Au coUisions at ^/s^ = 200 GeV have a 
dNch/drj similar to those of midcentral p+Ph collisions 
at the LHC, while the £2 values for d+Au collisions are 
about 50% larger than those calculated for the midcen- 
tral p-|-Pb collisions. The key observation is that the ratio 
^2/^2 is consistent between RHIC and the LHC, despite 
the factor of 25 difference in collision center of mass en- 
ergy. A continuation of this same trend is seen by also 
comparing to V2/£2 as measured in Au+Au [30l432| and 



Pb-hPb [33|, |3|| coflisions. 

In summary, a two-particle anisotropy at midra- 
pidity in the 5% most central d-l-Au collisions at 
y/s^ — 200 GeV is observed. The excess yield in cen- 
tral compared to peripheral events is well described by 
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FIG. 5: (color online) The eccentricity-scaled anisotropy, 
V2/S2, vs charged-particle multiplicity (dNch/drj) for p{d)+A 
collisions measured by PHENIX, ATLAS '^], and ALICE f|L 
Also shown are Au-f-Au data at ys„„ = 200 GeV [3CH33] 
and Pb-fPb data at y^s];;]; = 2.76 TeV [H,!!]. The V2 are 
for similar pr selections. Due to the lack of available multi- 
plicity data in p-|-Pb and d-f-Au collisions the dNch/drj values 
for those systems are calculated from HIJING [2g| . 



a quadrupole shape. The signal is qualitatively similar 
to that observed in long range correlations observed in 
p-|-Pb collisions at much higher energies, but with a sig- 
nificantly larger amplitude than that observed in 0-2% 
central p-|-Pb collisions at ATLAS. While our acceptance 
does not allow us to exclude the possibility of centrality 
dependent modifications to the jet correlations, the sub- 
traction of the peripheral jet like correlations has been 
checked both by varying the Ary cuts and exploiting the 
charge sign dependence of jet-induced correlations. The 
observed results are in agreement with a hydrodynamic 
calculation for c?-t-Au collisions at y/s„„ = 200 GeV. 

We find that scaling the results from RHIC and the 
LHC by the initial second order participant eccentricity 
brings the RHIC and LHC results to a common curve as 
a function of dNch / drj also shared by elliptic flow coeffi- 
cients from Au-|-Au and Pb-|-Pb collisions. This finding 
suggests that these phenomena are sensitive to the initial 
state geometry and that the same underlying mechanism 
is responsible in both p+Ph collisions at the LHC and 
d+AvL collisions at RHIC. It also suggests a relationship 
to the hydrodynamic understanding of V2 in heavy ion 
collisions. The observation of these correlations at both 
RHIC and the LHC provides important new information 
for understanding these phenomena. Models which seek 
to describe these features must be capable of also ex- 
plaining their persistence as the center of mass energy is 
varied by a factor of 25. 
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